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I. INTRODUCTION 



The absence of an interacting continuum limit or triviality of Abelian gauge theories in 
four dimensions such as QED is due to the screening of charged particles by their interactions 
with virtual fermion-antifermion pairs from the vacuum. Such charge screening is described 
by the vacuum polarization II (Fig. [I]). The QED vacuum is not a perfect insulator and 
can be considered as a medium of dipoles representing the fermion loops in the vacuum 
polarization. Within perturbation theory, the electromagnetic charge or gauge coupling is 
screened completely in the continuum limit (A — > oo, where A is the ultraviolet cutoff). 

This can be seen by considering the renormalization group (RG) transformation of QED 
H[, which relates the gauge coupling or fine-structure constant in the infrared (IR) region 
to the bare gauge coupling ao = oa in the ultraviolet (UV) region via 

a fl = H a [fi/A,a ], TZ a [///A, Qf ] = 1 - n ^ A a ^ ■ (1) 

In perturbation theory, the one-loop leading contribution to II is 

n[/VA,a ] = ^r^bi- (2) 

where N is the number of fermion flavors. This logarithmic screening effect is sufficient to 
cause the complete screening of charge in the continuum limit [0.0. 

From the RG point of view triviality is merely due to the absence of an UV stable fixed 
point or nontrivial root of the j3 function 



at \ - dTla K a °] 



dw 



w=l 



for the gauge coupling. Equation (0) gives rise to the j3 function (3 a {a) = 2Na 2 /3n, which 
only has an IR stable fixed point or trivial root: a = (/%(0) > 0). A nontrivial or 
interacting theory only arises whenever /3 a has a root which is an UV stable fixed point. In 
case of asymptotic free theories such as QCD the origin a = is an UV stable fixed point. 

In addition to its motivation only on the level of perturbation theory the above consid- 
eration is too naive. According to the RG methods of Wilson fH|, one should consider the 
RG flow in the space of all coupling constants (respecting certain symmetries) or at least 
for those coupling constants which can be classified as relevant close to a particular critical 
manifold in coupling constant space. For a possible nontrivial continuum limit of QED, 
critical and hence nonperturbative dynamical effects are required. 

Therefore, the discovery of dynamical chiral symmetry breaking (D^SB) in the strong 
coupling phase of QED (a > a c = 7r/3) |||| and the existence of an UV stable fixed point 
in the quenched-ladder approximation [pj-|iT|j sheds new light on the nonperturbative nature 
of QED and triviality. Lattice simulations of so-called noncompact quenched QED have 
confirmed the existence of a continuous chiral phase transition [O-O]. 



An important step was performed by Bardeen, Leung, and Love in Refs. [[T|,|T6|]. These 
authors realized that, in quenched-ladder QED, attractive four- fermion interactions de- 
scribed by a dimensionless four-fermion coupling go > have a so-called scaling dimension 4 
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instead of 6 at the critical gauge coupling «o = a c . Consequently, these operators mix with 
the gauge interaction which also has dimension 4 in four space-time dimensions; this means 
that QED is not a closed theory at the chiral phase transition. The scaling dimension of an 
operator is important in determining whether or not such an operator describes a long-range 
interaction. The model of Bardeen, Leung, and Love is referred to as the gauged Nambu- 
Jona-Lasinio (GNJL) model, and the mechanism of dynamical chiral symmetry breaking 
proposed by them should be regarded in the context of the RG methods of Wilson. 

Wilson pointed out that a nontrivial renormalizable model can only be formulated if the 
RG transformation exhibits UV stable fixed points. UV stable fixed points are "singular" 
points^ of the RG transformation at which the model becomes scale (conformal) invariant. 
Natural candidates for UV stable fixed points are critical points governing a continuous 
phase transition. Since, at a continuous phase transition, the correlation length £ is infinite 
and the model is scale invariant. 

The most crucial observation of Wilson is that in the RG transformation, i.e., the "coarse 
graining" process, new types of local interactions are generated and that the new interac- 
tions can be classified as either irrelevant or relevant interactions. Only relevant (including 
marginal) interactions are important in determining what kind of infrared dynamics char- 
acterized by a scale fi ~ l/£ emerges from the microscopic or bare model characterized by 
the cutoff A. The effect of irrelevant interactions can always be absorbed by adapting the 
coupling constants of relevant and marginal interactions. 

Especially, close to a continuous phase transition, the RG methods show that it is im- 
possible, a priori (without solving the equations of motion), to determine which interactions 
are relevant or irrelevant; particular interactions can acquire anomalous dimensions and in- 
teractions which are irrelevant in a certain region of coupling constant space might become 
relevant in another region. 

In this respect the GNJL model should be considered as the Wilsonian effective (or micro- 
scopic) action of QED taking into account the four-fermion interactions describing neutral 



scalar and pseudoscalar fermion-antifermion composites. It was shown in Refs. ||i7||18f that 
there is a critical line (curve) in the coupling constant plane (a , g ) of the GNJL model 
separating a chiral symmetric phase from the chiral broken phase. The critical line is given 
by 

g c {a ) = (1 + uj) 2 /A, u = y/l-AX , A = 3a /47r. (4) 

In the neighborhood of the critical line four-fermion interactions acquire sufficiently large 
anomalous dimensions to become relevant operators. The existence of nontrivial scaling 
behavior of the model near criticality implies that these scalar and pseudoscalar composites 
are relevant degrees of freedom at both short and long distances.0 

Lattice simulations of noncompact quenched QED with an induced four-fermion coupling 
were performed by the Illinois group in Refs. The Illinois group obtained a critical 

point (0.44a c , 0.76) in the (ao, go) plane, which fits nicely on the critical line Eq. (H). However 



1 UV stable fixed points are specific roots of the (5 functions. 

2 The composites are light states, as well in the symmetric phase as in the broken phase. 
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it was not possible at that time to investigate by means of lattice simulations the phase 
transition along the critical line. 

Later, the GNJL model was studied on the lattice in noncompact formulation using some 
mean-field approach for the fermions in Ref . |21j . They obtained a critical line qualitatively 



similar to the one following from the Schwinger-Dyson equation (SDE) approach Eq. 

In the intermediate region (0 < oq < a c ), the critical exponents describing the chiral 
phase transition satisfy nonmean- field hyperscaling relations which supports the view that, 
within the quenched-ladder mean-field approximation, the GNJL model has a nontrivial 
continuum limit [IT5,32[. Also in Ref. [21] nonmean field critical exponents were obtained. 



The physical implications and the consistency of the quenched-ladder results with many 
quenched lattice simulations, and with the nonperturbative RG techniques, support the view 
that the qualitative features of the approach might be realistic and describe properties of 
the full theory. Most likely this is due to the ladder approximation respecting the vector 
and axial Ward-Takahashi identities. 



In Refs. p3| , [24]| the validity of the ladder approximation was tested positively by includ- 
ing the effects of, e.g., crossed photon exchange graphs. In addition, the nonperturbative 
renormalizat ion-group methods of Refs. |25|j2^| provide a way to check the quenched-ladder 



approximation in the GNJL model by including the effect of crossed photon exchange graphs, 



and four-fermion interactions in the RG flow of coupling constants. In Ref. [26] the critical 
linef] and critical exponents in the full quenched GNJL model were obtained in a particular 
so-called local potential approximation, which incorporates besides crossed photon exchange 
graphs also four-fermion exchanges beyond the mean-field approach. Considering the small 
quantitative differences, qualitatively this study supports the reliability of the ladder ap- 
proximation. 

Attempts to include a logarithmic running of the coupling drastically changes the chi- 
ral phase transition and the critical line, see Refs. p7|-|30|. Moreover, it was shown in 
Refs. pp|-p5|] that the critical exponents are of the mean-field type (up to logarithmic viola- 
tions) leading to a trivial theory. 

Lattice simulations of noncompact full (unquenched) QED on the lattice (with flavors, 
N = 2 and N = 4) are controversial ||34)| . The Illinois group [0,0] (see also |3^|1^]) and 



the Zaragosa group |38| , |39| , find power-law scaling and nonmean-field critical exponents, sig- 
naling a possible nontrivial continuum limit for the strong-coupling broken phase, whereas 
p0| , |4T|j34]| obtain mean-field behavior (mean-field critical exponents with logarithmic correc- 
tions). Thus Gockeler et al. find a vanishing renormalized gauge coupling and a vanishing 
effective Yukawa coupling (defined by the Goldberger-Treiman relation), and they conclude 
that lattice QED is trivial, see for their most recent result Ref. [Q. 

In this paper we argue that the four-fermion interactions might play a crucial role in the 
phenomenon of charge screening. We show that, by considering the Abelian GNJL model, 
UV stable fixed points of the (3 function of the gauge coupling can be realized, provided the 
number of fermion flavors N exceeds some critical value. The existence of UV stable fixed 
points gives rise to a nontrivial theory. 



3 They obtained the critical line: g c (&o) = (l + ^) 2 /3, which differs by a factor 3/4 from the 
quenched-ladder result. 
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An important observation is that the nonmean-field values for the critical exponents sug- 
gest the existence of a nontrivial Yukawa interaction describing the interactions of the scalar 
and pseudoscalar composites with fermions. This also points out the inconsistency of the 
mean-field approximation (Hartree-Fock approximation) for the four-fermion interactions. 
Therefore, we go beyond the mean-field approximation by incorporating these composites 
(the a and it exchanges) in the 1/N expansion. 

The setup of this paper is the following. In Sec. [TT] we introduce the Abelian GNJL 
model with Ul(N) x Ur(N) symmetry. Furthermore, we sketch how we will search for the 
existence of an UV stable fixed point in the coupling constant plane (a ,g ). In Sec. |TJ 
we discuss, within the quenched-ladder mean-field approximation, the importance of hyper- 
scaling relations and how this is related to the existence of a nontrivial Yukawa interaction 
in the GNJL model. We discuss how to proceed beyond the quenched mean-field approxi- 
mation in Sec. |V[ The 1/N expansion is discussed in Sec. fV B\ In order to get some idea 
how the scalar and pseudoscalar composites contribute to vacuum polarization we illustrate 
such contributions on the level of perturbation theory for a gauge-Higgs- Yukawa model in 
Sec. |V[ Then, in Sec. [VT|, we argue how we can exploit the 1/N expansion and derive a 
computable expression for the (3 function, @ a , of the gauge coupling. Subsequently, (3 a is 
computed explicitly in Sec. |V11| . The existence of UV stable fixed points is addressed in 
Sec. |V111| . Finally, we present our conclusions in Sec. |TX|. 



II. THE ABELIAN GNJL MODEL 

We consider the GNJL model with U(l) gauge symmetry (the Abelian GNJL model) with 
N number of fermion flavors described by the following Lagrangian (see also Refs. P3| ,| 



^ a=0 



(5) 



where = <9 M — ie^A^ and where the flavor labels i, j run from 1 to N. The Lagrangian 
C\ is parametrized by the following three dimensionless bare coupling constants: 

/i = m /A = m A , a = e 2 /A7r = a A , g = G A 2 / An 2 = g A , (6) 

where A is the ultraviolet cutoff. We assume that the above set of dimensionless coupling 
constants comprises the entire set of relevant (including marginal) dimensionless coupling 
constants respecting particular chiral, vector, and gauge symmetries. 

If the bare mass mo is zero, C\ has a global Ul{N) x Ur(N) symmetry. The generators, 
r, of the U (N) Lie algebra have the following properties: 



N 2 -l 

r«t = T \ Tr t<*t? = S*, £ r«r« = 6 u S kj . (7) 



a=0 



The last identity is called the Fierz identity. 

The Abelian GNJL model described by Eq. (|5|) can be conveniently analyzed in terms 
of auxiliary or composite fields a a = —GQipT a if) and ir a = —Gq^t 01 (ij^ip describing scalar 
and pseudoscalar degrees of freedom. In this way C\ can be rewritten as 
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i N 2 -l _ i N 2 -l 

where mo has been set to zero. In this formulation the four-fermion interactions are described 
by the interactions of the auxiliary fields with the fermion fields. Then the connected two- 
point Green functions of the ir a fields describe N 2 Nambu-Goldstone bosons (tt bosons) and 
the connected two-point Green functions of the a a fields describe iV 2 "Higgs" bosons (a 
bosons) . 

Although the Lagrangian (||) has very interesting properties, surprisingly this particular 
class of GNJL models with Ul(N) x Ur(N) has not received much attention. One such 
property is that Eq. ([5]) comprises the largest set of relevant chiral invariant four-fermion 
operators for a model of N fermions. Another independent set of vectorlike chiral invariant 
four-fermion interactions such as (^7^) 2 do not acquire large anomalous dimensions and 



remain irrelevant near the critical point (line). This has been shown in Refs. p5| , |26|| , where 
the RG flow of scalarlike chiral invariant and vectorlike chiral invariant four-fermion inter- 
actions has been considered. Moreover, the specific form of the chiral symmetry, where the 
number of scalars equals the number of pseudoscalars, turns out to have useful implications 



in the context oil/N expansions as we will discuss later in Sec. |TV. 
In case of £2 there are four renormalization constants: 

Z 2 1/2 (AVA)^ (A0 (x) = if>(x), zl'\A'/A)A^ l) (x) = A»{x), (9) 
Z a (A'/A)a [AI) (x) = <r(x), Z„(A' /A)tt {AI) (x) = n(x), (10) 

where A'/A < 1, and the fields if), A^, tt, a are the bare fields defined at the UV cutoff A. 

Near the critical point a new scale is generated: the correlation length £. In case of a 
second-order type of phase transition, the correlation length exists in both phases of the 
system. In the broken phase the inverse correlation length is real and can be considered as 
a physical mass of particles, for instance the mass of the scalar bound state (the a boson) 
ma, or the mass of fermion md yn - In the symmetric phase the fermion is massless, and the 
scalar and pseudoscalar composites are unstable states characterized by a complex mass 
pole in their respective propagators describing the mass and the width of the Breit-Wigner 



type resonance, see Refs. [^5 -|^|. The absolute value of the complex mass pole \m a \ can be 
considered as the inverse correlation length, i.e., \m a \ ~ 

The RG transformation dictates the flow of the dimensionless bare couplings as function 
of the UV cutoff A. Typically the bare relevant and marginal couplings (e.g., /i , a , g ) have 
to be fine-tuned sufficiently close to the critical point in order for scaling behavior to set in, 
so that the physics in the infrared can be related to experimental data. Scaling behavior is 
obtained when there is a large scale hierarchy between the infrared length scale £ and the 
ultraviolet length scale a — 1/A, i.e., £ ^> a. 

The fine-tuning depends on the eigenvalues of the RG transformation of the couplings 
close to the critical point and hence on the critical exponents. These critical exponents can 
be derived from the f3 functions for the coupling constants /x , ot Q , and g ; 

A-^ = Pp(fio, 010, go), A^ = f3 a (no, a ,g ), A^ = /3 g (/i Q , a ,g )- (11) 



The crucial step is determine the fixed points (//*, a*, of the RG equations fllTD 



i.e. 
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/^(/x*, a*, g*) = 0, P a {fu, a*, #*) = 0, /? 9 (/x*, a*, gf*) = 0, 



(12) 



since the nature of the fixed point determines whether a nontrivial continuum limit can be 
realized or not. For a nontrivial continuum limit (/!*,«*,#*) should be a UV stable fixed 
point. 

The RG equations follow from the regularized SDE's of the generating functional de- 
scribed by the Lagrangian (|]).f] The UV stable fixed point for (3^ is /x* = 0, hence we can 
write 

A^^-(l + 7m )/xo, (13) 

where 7 m is the anomalous dimension of the mass operator ipip evaluated at the fixed point 
(/x*, a*, g+). In the case of the quenched GNJL models 1 < 7 m < 2. Thus the dimensionless 
bare mass /xo is a relevant coupling requiring fine-tuning. 

After setting /x = /x* = 0, the problem reduces to the determination of the UV stable 
fixed points in the coupling constant plane (a ,g ): i.e., 



Pa{<x*,g*) = 0, 



dp g (a,g) 

dg 
d/3 a (a,g) 



= 0, P' g (a*,g*) = 9 " — <0, (14) 

dg (a, S )=K, S *) 

< 0. (15) 



(9a 



(a,s)=(a*,ff*) 



The quenched-ladder approximation simplifies the solutions of Eqs. (|T4| ) and ( |T5"D con- 
siderably since the quenched hypothesis explicitly sets /3 a = for all «o by omiting fermion 
loops. It was shown in Ref. ||8| that, in the symmetric phase (go < g c ), 

/3 g (ao, go) = -2u—(g - g c ), (16) 

gc 

with uj and g c given in Eq. (Q). In the next section, Eq. ( p!6|) will be derived. Clearly, in this 
way, the UV stable fixed point of j3 g is the critical line; 

g* = g c (a ), (3 g (a , g c (a )) = 0. (17) 

Now Eq. ( |15"D should be reconsidered. We will analyze Eq. (|15[) beyond the quenched ap- 
proximation, and try to solve 

Pa(ao, g c (a )) = 0. (18) 

In Sec. |VT| an explicit expression for f3 a will be derived by assuming that go is at its critical 
value g c and that it has an UV stable fixed point so that in the neighborhood of this point 



A derivation of the set of (full) SDE's for the two and three-point functions is given in Chap. 2 
of @. 
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III. HYPERSCALING IN THE QUENCHED-LADDER MEAN-FIELD 

APPROXIMATION 



In analogy with statistical mechanics, the continuous chiral phase transition can be 
classified in terms of critical exponents which describe the scaling of various macroscopic 
quantities (e.g., the chiral condensate, correlation length, effective potential, chiral suscep- 
tibility) close to or at the critical point. It is considered a strong indication of the existence 
of a nontrivial continuum limit (A — > oo), if so-called hyperscaling relations between these 



various critical exponents are satisfied, see Refs. flT||22| . |49| - 

Because the a boson propagator A5 is the connected correlation function of the field 
a describing correlations parallel to direction of symmetry breaking (i.e., parallel to the 
direction of long-range ordering), the absolute value of the mass, m a , of the a boson, which 
is given by As, is the natural candidate for the inverse correlation length. 

In the quenched- ladder approximation treating four-fermion interaction in a mean-field 
approximation, the critical exponents are 

r 2 + ^ 2-uj 1 

d =o > = -^ — ' U= 7T> 7 = 1, 19 

2 — uj 2uj zoo 

and satisfy the hyperscaling relations 

y = 0(6-1), 4i/ = 2/3 + 7 . (20) 

Other hyperscaling relations involving the critical exponent a describing the scaling of the 
effective potential are satisfied too Furthermore, it was argued in jHJ that the 



validity of the quenched-ladder mean-field approximation relies on the verification that the 
critical exponent 7 = 1. The interpretation of 7 = 1 is the factorization 77^)2 = 2r]r^y 
A renormalization of the chiral condensate simultaneously renormalizes the propagators As 



and Ap. Indeed the lattice computations of the critical exponent 7 reported in [|T9|-|21 
showed strong evidence for 7 = 1. 

The anomalous dimension 77 describes the scaling of the connected two-point Green 



function As at the critical point. In Ref. |47] the scalar propagator As (see Fig. ||) and the 



scalar Yukawa vertex Ts (see Fig. |3|) have been computed in the symmetric phase in the 
quenched-ladder approximation by means of a so-called two-channel approximation. 
In the symmetric phase the Yukawa vertex has the following form: 

T s (p + q,p) = l[F 1 (p + q,p) + (qp-pq)F 2 (p + q,p)]. (21) 

These vertex functions F% and F2 can be expanded in terms of Chebyshev polynomials of 
the second kind, e.g., 

00 P • Q 
Fi(p + q,p) = V f n (p 2 ,q 2 )U n (cosa), cosa = . (22) 

^0 PI 



The two-channel approximation of Ref. [17] now holds in that the Yukawa vertex is approx- 
imated by the angular average of the vertex function F\ in the following way: 

T s (p + q,p) « 1 J < ^F l (p + q,p) = lf (p 2 ,q 2 ), (23) 



S 



where 



W 2 , q 2 ) = F m {p\ q 2 )d{q 2 - p 2 ) + F uv (p 2 , q 2 )9(p 2 - q 2 



(24) 



The functions Fir, and Fuv are, respectively, referred to as the IR channel (infrared), and 
the UV channel (ultraviolet). The specific choice of Chebyshev expansion and consequently, 
the choice of zeroth-order coefficient /o, is convenient, since the infrared limit (q 2 ^> p 2 ), 
and ultraviolet limit (q 2 <C p 2 ) of Y s are both described by f . The scaling form for As(q) 
is well described by these two functions F IR and Fuv- 

The zeroth-order Chebyshev expansion or two-channel approximation gives second-order 
differential equations for the lowest order harmonic with appropriate IR and UV boundary 
conditions. These differential equations are exactly solvable and the solutions are expressed 
in terms of a Bessel function of the first kind for Fir and in terms of modified Bessel functions 



for Fuv- The solutions are given in pTFfl . 

With the "asymptotic solutions" for Ts given in terms of Fir and Fuv, an analytic 
expression for the a boson can be obtained from 



4^2 Y Q 



F uv (A 2 ,g 2 



As X (?) = - 



1 



n s (g 2 



(25) 



which has been argued in [|7| to be correct up to leading and next-to- leading order in q 2 /A 2 
and is valid along the entire critical curve g c . 

Also one can derive that the solutions for Ag and are consistent with hyperscaling for 
< a < a c . This is intimately related to the fact that the renormalization of the auxiliary 
fields a and tt, Eq. (|10D, simultaneously renormalizes the Yukawa vertex and the scalar 
propagator. 

Near go = g c (with \m a \ 2 ) q 2 A 2 ) the scalar propagator has the scaling form 
44,49-[5TJ (in Euclidean formulation): 



m ? ) = ^ (f 



l-ri/2 



-Fa(KI7<7 2 ) 
where i] is the anomalous dimension 

77 = 2(1- 



•FaOe) 



4tt 2 



1 



B(uo) l + x u 



and where 



B(u) 



16a; 7(-w) T(2 - w) 



1-lu 2 ) 2 7 H r(2 + w) 

2AT)/I 



2A 



2 A, 




(26) 
(27) 

(28) 
(29) 



With Eq. (|T9D, the anomalous dimension rj, Eq. (p7|), satisfies the hyperscaling relation 

7 = 1/(2-7;). (30) 



It has been shown in [E5 



that As has a complex pole on a second Riemann sheet in 
(Minkowsky) momentum space. The mass m a is the complex pole of As and the absolute 
value \m a \ scales according to \m a \ ~ (Ag ) u . More precisely, in [^5] , ^7| it is derived that 
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ImJ oc A 



g c g B(Lu) 



l/2w 



Ag = go- g c - 



(31) 



The absolute value of m a is taken to be the physical, macroscopic or infrared mass scale, 
which by definition should be independent of A. Using Eq. (0), we can derive the (3 function 
of go by assuming that 



= A 



d\m a \ 
dA 



(3 g (aio,go) = ~^—(go ~ 9c) 



(32) 



which is equivalent to the j3 function given in Ref. ||8| and Sec. ||. Hence the critical curve 
go — gc is an UV stable fixed point fi g (ao,g c ) = (/3' (ao,g c ) < 0) of the RG flow. 

In accordance with Ref. j|7[], the scaling form for the Yukawa vertex can be written as 



T s (p + q,p) » 1 ( ^ ) [Jth (p 2 /q 2 ) 6(q 2 - p 2 ) + F vy (q 2 /p 2 ) 6(p 2 - q 2 



2\ V/4 



(33) 



where, for p , q <A , 



Fm(p\q 2 ' 



'A 



2\ 



'A 



2\ 



(34) 



and 



Fw(q 2 /p 2 ) , . ... 

V / 7(0;) (1 + u;) 



2sin^7r 2 r(l -w) fX Y /2 fq 2 V /2 



7T 7(to>) (1 + W) 

2 



'2A p 2 



(35) 



2 r(l-w) A \" /2 /V\ 1/2 



7(w)Z- 



/2A g 2 



-7(-o;)4 



'2W 

p2 



(36) 



From the above scaling form for As and it is clear that four-fermion scattering amplitudes 
such as 



r<?(pi + q,p 1 )A s (q)T s (p2,P2 + q) oc — , p 1? p 2 < g <A 



(37) 



are independent of A and express the long-range natureQ of Yukawa forces. For dimensions 



2 < d < 4, this was pointed out in Ref. ||46|| . The long-range Yukawa forces and their 
nontrivial contributions to scattering amplitudes in the infrared are a direct consequence of 
hyperscaling and thus powerlike renormalizability. 

The consensus is that in four dimensions due to the logarithmic corrections, the hy- 
perscaling relations are violated for the pure NJL model and A0 4 theory, and we have the 
following inequalities: 



5 The correlation length is large, l/|m CT | S> 1/A. 
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4z/>2/5 + 7 , j > (3(5-1), (38) 



see Ref. []52[ for an extensive discussion regarding this issue. This violation of hyperscaling is 
believed to be a sign of triviality meaning that the effective Yukawa coupling (which couples 
Goldstone bosons to the fermions) vanishes in the continuum limit. The continuum limit 
is noninteracting, hence trivial. This can be seen in the following way. Assuming that in 
the low-energy region the correlation length £ is the only relevant length scale, we define 
an effective Yukawa coupling gy by the zero-momentum limit of the four-fermion scattering 
amplitude (two fermions exchanging a scalar bound state) in the D^SB phase 



2 



% ~ e g 2 Y ~ r s (o, o)A 5 (o)r 5 (o, o), (39) 

where r,s(0,0) and A s (0) are given by the chiral susceptibility relations 

where So = S(0) is the fermion mass in the broken phase. By making use of the scaling 
laws JI9] , P^| and that £ ~ 1/S ~ \jm a , (er) ~ (ipip), we can derive that 



g y ~ £(2/3+7-4,)/^ (41) 



This expression is related to the definition of gn ~ gy given in Ref. and it is clear that 
the scaling inequalities ( j38|) imply that g\ — > when £ — > 00. Thus, only if the hyperscaling 
relations (^) are satisfied, a nonzero gy might be realized in the continuum limit (£ — > 00), 
thereby giving rise to a nontrivial interacting theory. 



IV. BEYOND THE QUENCHED-LADDER MEAN-FIELD APPROXIMATION 

The quenched approximation is analogous to the assumption that the full photon prop- 
agator 

o„M = (-9,^ + q ~f) AM - ^ ^Y^y (42) 

can be approximated by the bare or canonical propagator A(q) = 1/q 2 (for large momenta), 

DM = (-9^ + qj %) 4, (43) 



in the Landau gauge (a = 0). The quenched approximation is only consistent when the 
vacuum polarization is finite in the continuum limit, i.e., the logarithmic running of the 
coupling is absent. This is the case at an UV stable fixed point of the (3 function, Eq. (|15|), 
of «o- The assumption that such a critical fixed point exists, and that it lies somewhere on 
the critical curve (H) is the starting point for many studies of dynamical chiral symmetry 



breaking in context of the GNJL model. In fact, the quenched hypothesis (H3I) is only 
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consistent when the bare gauge coupling a is near the fixed point of the theory, (3 a ~ 0. 



We discuss this issue in more detail in Sec. VI 



In many approximations of the GNJL model, the four-fermion interactions are treated in 
a mean-field approach known as the Hartree-Fock approximation. In mean-field approxima- 
tions the composite operators such as ipip are replaced by their vacuum expectation values 
(-0-0 — ► (ipip)) and fluctuations about that value are ignored. Thus quantum corrections 
corresponding to four-fermion interactions are neglected beyond tree level. 

As long as four-fermion interactions are irrelevant the mean-field approach for these 
operators is justified. However, in Refs. [0,^] it is concluded that, within the quenched- 
ladder mean-field approximation to the GNJL model, the hyperscaling equations for the 
critical exponents are satisfied, implying that the four-fermion operators become relevant due 
to the appearance of large anomalous dimensions. In other words, the mean-field approach 
yields non-mean-field exponents, thereby being inconsistent (e.g., see Refs. [[H|,|5(J). As was 
discussed in Sec. PH the hyperscaling relations imply the existence of a nontrivial Yukawa 
interaction describing the interaction between fermions and a and tt composites in the GNJL 
model at both short and long distances. 

In order to go beyond the mean-field approach, we propose the following. First, we point 
out the usefulness of skeleton expansions. Second, we make use of the specific form of the 
chiral symmetry and adopt the 1/N expansion (with iV the number of fermion flavors). 



A. The skeleton expansion 

The non-mean-field values of the critical exponents imply that one cannot neglect (as is 
done in mean-field approximations) the full connected Green functions corresponding to the 
composites (or at least the leading or asymptotic parts of these functions) in the SDE's. 

On the level of the Bethe-Salpeter (BS) fermion-antifermion scattering kernels the a and 
tt composites can be incorporated in a RG invariant manner by making use of the skeleton 
expansion, e.g., see Ref. [[53]). Analogous to QED kernels we define the one-boson irreducible 



kernel and the two-fermion one-boson irreducible BS kernel K^ 2 \ where these kernels 
now also include the a and tt composites. For both type of kernels a skeleton expansion 
exists. The integral equation between and is known as the Bethe-Salpeter equation. 

The skeleton expansion is a series in topologically distinct Feynman diagrams with all 
vertices and propagators fully dressed. The skeleton expansion is a special way of resumming 
the entire set of Feynman diagrams in a consistent manner, i.e., without double counting. 
The lowest order terms ("lowest" in terms of loops) of the skeleton expansion for is 
illustrated in Fig. The blobs with the letter "B" in the full vertices and propagators 
represent photons, and composite a and tt exchanges. 

Each term in the skeleton expansion of the BS kernel is RG invariant, up to fermion 
wave function factors, i. e. the expansion is independent of the renormalization factors Z 3 and 
Z = Z a = Z w (see Eqs. (||) and (|10D ) of, respectively, the gauge field and the composite fields 
a and tt. The two Z~ l factors with anomalous dimensions of each Yukawa vertex cancel with 
the Z 2 factors of the a and tt propagators, leading to cutoff independent fermion-antifermion 
scattering amplitudes, e.g., see Eq. fl3"T|). 
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B. The Ul(N) X Ur(N) chiral symmetry and the 1/N expansion 



As was mentioned in Sec. IOL the Abelian GNJL model, with iV number of fermion 
flavors, is taken to be invariant under global Ul{N) x Ur(N) chiral transformations, so that 
both the scalar and pseudoscalar four-fermion interactions are in the adjoint representation, 
and, consequently, the number of scalar composites (iV 2 ) equals the number of pseudoscalar 
composites {N 2 ). In this way, when N is large we can use the 1/N expansion introduced 
by 't Hooft |54j]. This provides us with a scheme to incorporate four-fermion interactions 
beyond the mean-field approach. The 1/N expansion states that the planar (i.e., ladder) 
diagrams, with fermions at the edges, describe the leading or dominant contributions to 
Green functions. 

The interesting feature of such a 1/N expansion is that Feynman diagrams can be classi- 
fied in terms of two-dimensional surfaces with specific topology. Diagrams with other (than 
planar) topological structures are suppressed by at least factors of 1/N, and in the limit of 
large N, their contribution can be neglected with respect to planar graphs. One important 
rule is to draw Feynman graphs with fermion loops forming the boundary of the graph (if 
possible). In this way, vertex corrections are not necessarily classified as being planar. 

In the context of the 't Hooft's 1/N expansion, we should consider internal or virtual a 
and it exchanges analogous to the gluon exchanges with the important difference that due to 
the chiral symmetry we have two types of particles both being in the adjoint representation 
(iV 2 scalars and N 2 pseudoscalars) . Then, by keeping track of the flavor indices within a 
particular Feynman diagram, we can count factors of 1/N. Each fermion carries a flavor 
index (i), which runs from 1 to N. A virtual a, tt exchange, being associated with two 
Yukawa vertices, carries two flavor indices. Therefore, as a result of the Fierz identity (0), 
each virtual a, n exchange gives rise to a pair of Kronecker 5 functions connecting the flavor 
indices of the scattered fermions. In the context of flavor indices, either a a or it boson can 
be considered as a propagating fermion- ant if ermion pair carrying double flavor indices. 

Whenever a trace over a flavor Kronecker 5 function enters into the expression for a 
particular Feynman diagram, we speak of an index loop. An index loop is easily identified 
by using the double-line representation of 't Hooft. A fermion propagator is represented by a 
single index line (i.e., fermion line), whereas each internal scalar, respectively, pseudoscalar 
propagator is represented by a double index line. Consequently, whenever, after drawing a 
particular Feynman diagram, an index line closes, it forms an index loop giving rise to a 
factor N = Tr 5. 

The topology of the Feynman diagram determines the factors of N. The vacuum polar- 
ization has the topology of a sphere with a single hole (i.e., a disk), where the fermion-loop 
forms the boundary (i.e., hole) of the graph. It can be shown straightforwardly, that planar 
diagrams in the vacuum polarization with n exchanges of ex's and 7r's are associated with 
a factor iV n+1 gy n , where gy is an "effective" Yukawa coupling describing the interaction 
of scalars and pseudoscalar with fermions. For the time being we leave unspecified such a 
coupling. 

In absence of bare mass, the Ul(N) x Ur(N) symmetry allows us to write each full 
Yukawa vertex, photon-fermion vertex, fermion propagator, and o, ir boson propagator as 

T a sab (k,p) = T«T Sab (k,p), A^(g) = A s (q), (44) 
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T a pab (k,p) = r«T Pab (k,p), A { °\q) = A P (q), (45) 



T%(k,p) = 5 l3 T^ b (k,p), S®(p) = S(p), (46) 



with a, b spinor indices, i, j flavor indices, and a the [/(iV)-generator index (see Refs. HJT7 



for the definitions of the proper vertices and connected two-point Green functions). So that 



N 2 -l 



a=0 



£ T a sab (k + q, k)A { « ) (q)T% d (p,p + q) = 5 u 5 kj T Sab (k + q, k)A s (q)T Scd (p,p + q), (47) 



because of the Fierz identity (0). Then the first term of the skeleton expansion for is 
the following single boson exchange term: 

{-iel)Kf b } cd (k,p,p + q) = 5 ilh 5 i2j2 (-i)Y S cb(p + q,p)iA s (q)(-i)T S ad(k,k + q) 

+ tiiijAihi-i)^ 'pcb(p + q,p)iAp(q)(-i)T Pad (k, k + q) 

+ Si^6i da (-ieo)r^(p + q,p)iD}^(q)(^eo)I^(k,k + q). (48) 

As a result of the chiral symmetry the contributions of four-fermion interactions, which 
are represented by a and it exchanges, exhibit two distinct features depending on whether 
they are incorporated in SDE's describing quantities connected with so-called zero-spin 
structures^ (e.g., the dynamical mass S, the Yukawa vertices Ts, Tp, and the a and n prop- 
agators As, Ap), or whether the exchanges are included in SDE's describing nonzero-spin 
structures (anti-commuting with 75) (e.g., the vacuum polarization II, the photon-fermion 
vertex T^, and the fermion wave function Z = Z 2 ). Henceforth, we refer to (non)zero-spin 
functions, and their equations as (non)zero-spin channels. 

The chiral symmetry gives rise to the following properties. 

1. In spin-zero-channels, the contribution of planar diagrams (i.e., planar in a and 7r 
exchanges) vanishes due to the fact that the exchange of a a has an opposite sign 
with respect to a tt exchange. Why? Let us consider a planar contribution to the 
scalar vacuum polarization which contains (amongst others) a 7r exchange. Both 75 
matrices corresponding to this particular planar tt exchange can be eliminated from 
the fermion trace of the scalar vacuum polarization by moving them to the right-hand 
side of the trace. For planar diagrams such a process involves the interchange of the 
75 matrix with an even number of fermion propagators, and an arbitrary number of 
Yukawa vertices. Since the Yukawa vertices commute with the 75 matrix, and 75 anti- 
commutes with the fermion propagatoij] S, the process of moving the 75 to the right 
does not introduce an overall minus sign. Now using that (275) (275) = — (1)(1), we 
see that the diagram containing a specific planar tt exchange is identical to minus the 
same diagram with the tt exchange replaced by a a exchange. Since each diagram 
containing a tt exchange has a scalar counter part (i.e., an analogous diagram with a 



6 Such structures are characterized by spinor matrices which commute with the 75 matrix. 
7 In the symmetric phase 75 S = —S^. 
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a instead of a 7r exchange), the sum of all planar diagrams, with a particular number 
of exchanges, vanishes. 



2. In nonzero-spin channels (think of II, T M , etc.) containing vertices which anti-commute 
with the 75 matrix, the situation is different: planar a and 7r exchanges contribute 
with identical sign. Let us now consider a planar contribution to the (photon) vacuum 
polarization containing a it exchange. If we again move the 75 matrices to the right- 
hand side of the trace, we get an overall minus sign due to the anti-commutation of 75 
with 7 M . This means that any planar diagram in the vacuum polarization containing 
a 7r exchange is identical to the same diagram with the 7r exchange replaced by a o 
exchange. 

The properties described above are, strictly speaking, only valid in the symmetric (massless) 
phase, where the (2N 2 ) a and tt bosons are degenerate. However, in the broken phase, 
the properties are valid whenever momenta larger than the dynamical mass £ or m CT are 
considered, because then the degeneracy emerges too. 

These properties also provide us with a general argument why the mean-field approach 
for four-fermion operators for Green functions corresponding to spin- zero channels (e.g., 
Ts and As) is reliable. For such channels planar contributions vanish and the next non- 
vanishing contributions (such as contributions containing crossed a and it exchanges) are 
proportional to l/N, thus small for large N. This implies that quantities such as the critical 
curve, dynamical mass, anomalous dimensions etc., are nearly independent of N, and are 
described rather well by the mean-field approach. To the contrary, the cancellation of scalars 
against pseudoscalars degrees of freedom does not occur in the vacuum polarization II which 
is a nonzero-spin channel. 



The inclusion of relevant four-fermion interactions beyond the mean-field approach re- 
quires a reinvestigation of the SDE for the fermion wave function Z = Z2 (Eq. @), 



with S the fermion propagator and S the dynamical mass. In QED in the quenched-ladder 
approximation, the fermion wave function has a gauge dependent anomalous dimension. In 
the Landau gauge, this anomalous dimension vanishes and Z — 1. 

We conjecture that the inclusion of relevant four-fermion interactions does not introduce 
an anomalous dimension for the fermion wave function other than already introduced by the 
gauge interactions. Thus, in the Landau gauge, the wave function Z is finite though it might 
deviate from unity. The main argument in support of the conjecture stated above is that 
only one full Yukawa vertex appears in the self-energy part, which means that anomalous 
dimensions of four-fermion interactions are not canceled. Only two fully dressed Yukawa 
vertices and a fully dressed scalar composite are RG invariant (anomalous dimensions can- 
cel). Consequently a remnant power of the cutoff (related to anomalous dimension of a 
Yukawa vertex) lowers the degree of divergence of the self-energy part from a logarithmic 
divergence to a finite integral. Therefore, throughout this paper we assume Z — 1. 



C. The fermion wave function 



S(p) 




(49) 



p — S(p 2 ) ' 
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A nice feature of the assumption that Z = 1 is that with the gauge interaction treated 
in the quenched-ladder approach the chiral and vector Ward-Takahashi identities (WTI's) 
are preserved, since in channels with spin-zero the planar a and 7r exchanges cancel each 
other. 



V. SCALARS, PSEUDOSCALARS, AND CHARGE SCREENING 

Since, in the GNJL model, the scalars and pseudoscalars are neutral states which there- 
fore do not couple to the photon field, their contribution to the vacuum polarization is 
described indirectly in terms of photon-fermion vertex corrections, and fermion self-energy 
corrections. Hence, in order to gain some intuition for the role of scalar degrees of freedom 
on the mechanism of charge screening, we analyze the two-loop contribution arising from a 
and 7T exchanges to the vacuum polarization. 

Let us consider an Abelian gauge-Higgs- Yukawa type of interaction described by the 
LagrangianQ 



e V>7 M ^V/> _ 9y${(t + i^)^ ~ V(a, it), (50) 



where the potential V contains, e.g., mass terms, and a a 4 type of interaction (i.e. a quartic 
scalar interaction). For simplicity, we ignore the effect of the potential V. In Appendix [A|, 
the two- loop contribution to II has been computed for the special case of N = 1, see also 
Fig. |5|. If the scalar and pseudoscalar fields in Eq. (|50|) are both in the adjoint representation 
of U(N), the result, for arbitrary N, reads 

n [q/A, a , X Y ] w — ^ f - + ^ - In - + {a /ir)O{l), (51) 

with Ay = gy/^Tr. The (3 function corresponding to such a vacuum polarization (Fig. H) is 

The interesting result of this computation is difference in sign between terms corresponding 
to photon exchanges, and terms corresponding to (pseudo)scalar exchanges. From this point 
of view, the fundamental scalars and pseudoscalars in a gauge-Higgs- Yukawa system tend to 
decrease charge screening. Furthermore, we might be tempted to conclude that a nontrivial 
root of Eq. (^) could be realized whenever NXy/2tc ~ 2/3. However, the complete situation 
is more involved. The RG equation for, e.g., Ay should be considered too, i.e., we should 
compute the (3 functions of Ay, and of any quartic scalar coupling. If and only if a nontrivial 
(nonzero) UV stable fixed point for Ay exists, the realization of a zero of Eq. ( |52D becomes a 



8 For a discussion of the renormalizability of non-Abelian gauge-Higgs- Yukawa models and non- 
Abelian GNJL models we refer to Refs. Ip5fl, and references therein. 
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realistic option. In other words, such a scenario is only possible if the Yukawa interaction Ay 
is nontrivial. The discussion in Sec. Ill implies that in order to obtain a nontrivial Yukawa 
coupling, the hyperscaling laws should be obeyed. The idea is that, nonperturbatively, close 
to the critical curve in the GNJL model, the scalar and pseudoscalar Yukawa interactions 
are nontrivial, and kinetic terms for the scalar and pseudoscalar composites are effectively 
induced via the appearance of a large anomalous dimension. 



VI. THE VACUUM POLARIZATION IN THE 1/JV EXPANSION 



The purpose of the present paper is to investigate the existence of an UV stable fixed 
point of the gauge coupling, so that 



0. 



Va = 0*) > 0, 



(53) 



with r] a a critical index characterizing the RG flow in the neighborhood of («*,#*). With a 
close to a+ the (3 function linearizes as 



P a {<Xi9*) ~ V a (ot* - a). 



(54) 



The (3 function for the gauge coupling follows from the RG transformation relating the 
charge a x at scale x to the charge a y at scale y via 



a, 



Oy 

/y,<Xy,g v ] 



II [l,a x ,g x ] = 0, 



(55) 



where IT(x 2 ) = II [x/y,a y ,g y ] is the (Euclidean) vacuum polarization. The RG transforma- 
tion should satisfy the RG semigroup property (with x < y < z) 

TZ a [x/y,a y ,g y ] = TZ a [x/z,a z ,g z ] = TZ a [x/y,K a [y/z, a z , g z ] , TZ g [y/z,a z ,g z ]} , (56) 
Tla[x/y,at*,g*\ =fc a lx/z,a±,g±], (57) 

where TZ g is the RG transformation for g, satisfying analogous equations. Then the (3 
function for a is defined as 



Pa is^-x i gx) 



dTZ a [w,a x ,g x 



dhiw 



(58) 



w=l 



The RG semigroup property ([56]) gives rise to the well known differential RG equation 

dot 

Pa(a y ,g y ) =y~j JL , ( 59 ) 



whose solution is 



a 3 



n=0 



7V. 



Pa {ptyi gy) 



d 
da,, 



a. 



(60) 



To obtain a nontrivial theory in the IR the existence of an UV stable fixed point is required. 
Close to the UV stable fixed point (P a — P g — 0), with g y — g+, we have that 
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Pa(ot y ,g*) 



d 



da 



y J 



a y « /? a (a y , (/%(<**, </*)) 



n-l 



thus 



'/a 



(61) 



(62) 



which is the solution of Eq. (|59|) using Eq. fl54]). The above expressions are only valid if both 
a x and a y are in the neighborhood of the UV stable fixed point a±, therefore we have the 
fine-tuning condition: 



(a y - a*) fy\ Vo 



x 



< 1, 



(63) 



which in case of y 3> x implies that a y is tuned sufficiently close to a*. 

In case of the Abelian GNJL model the gauge coupling a q in the infrared (IR) region is 
related to the bare gauge coupling a via the RG transformation a q = 1Z a [q/A, a , g ], where 
the four-fermion interactions contribute to the vacuum polarization IT (Fig. [I]) through the 
full photon fermion vertex, see Fig. || Assuming that the hyperscaling laws are satisfied, we 
can write the vacuum polarization for A ^> q > \m a \ (i.e., near g c {ceo)) as 



U[q/A,a ,g ] = /(i)(a ,. 



In — + /(2) (ao, go) f ln "~j +/(3)(«o,») ( ln "^) 



+ . . . 



(64) 



where each factor of InA/g corresponds to a single fermion loop with two outgoing photon 
lines. Thus the function /m represents the contribution of all diagrams to II in which the 
internal photon propagators are replaced by the bare or canonical form. 

In order for IT to give rise to a RG transformation satisfying Eqs. ( (55]) and (|56D, it can 
be derived from Eq. ( pD| ) that the functions , /( 2 ) , and /( 3 ) should be related to (3 a in the 
following way: 



®f{i)(ac,g) = Pa(a,g), 

f(l) ~ f(2) = 



2a 



/(l) - 2/(1) f(2) + /(3) = 



oa 



1 2 



Pa 



(65) 
(66) 

(67) 



These identities are nontrivial and require a high degree of self-consistency of the theory in 
the form of Ward identities.^ Since, within our approximation scheme, the Ward identities 
are respected, we assume that Eqs. (|65|)-(|67D are satisfied. Equation flS5[) now relates the (5 
function to the function f(\y Clearly the f3 function has a Gaussian or trivial fixed point at 
a-o = 0. 

From this we also have that with a — > a+ at g Q = g+ (see Eqs. (|54D , (162]) , and (|63])): 



9 The proof of Eqs. (|65|)-(|67|) can be performed order by order within perturbation theory. 
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II [g/A,a ,#* 



Hence, A(g) of Eq. (|^) is 




A( ff )«l(l-£^ 




<C 1. 



(68) 



-1 + 



.1. 



I la ■ 



(69) 



The second term on the right-hand side of Eq. fl69|) will only contribute via internal photon 
propagators to the functions /( 2 ), /p), etc. and not to the function /m. In order to find 
an UV stable fixed point of (3 a , we only have to compute /m and, therefore, we neglect 
all corrections to A(q) other than canonical in internal photon propagators. Because the 
contribution of II is neglected, this procedure is identical to quenching the internal photon 
propagators, although it is not quenched in the sense of taking iV — > 0! 

Now in correspondence with Eq. (|18|) we search for UV stable roots of Eq. (|65|), i.e., 
/(i)(«0) 9c(oio)) = 0. In case of pure QED (without four-fermion coupling) the function /(i) 
has been studied thoroughly by Johnson et al. in Refs. |56|-|58| and by Adler |59|] in the 
context of massless QED. In Ref. fl56 |, an expression for /m is obtained in term of the BS 



kernel (Sec. [TV]) as the single unknown Green function. 

We mention that, although the strong belief of Johnson et al. in the possible existence 
of finite QED seems poorly motivated from the point of view of Wilson's RG methods [f|.[!3 
their methods and techniques are sound and directly applicable to the GNJL model. 

In Appendix [B], we expose a brief derivation of the result of Ref. |56| and point out 



the applicability to the Abelian GNJL model. Then, by taking into account also relevant 
four-fermion interactions at go = g c (oto) via the BS kernel, we can derive from Eqs. (|B1|) and 
(m) that 



Na.% 

7i 



+ 



(70) 



where the functions (pi, <p2, and 03 are defined as follows: 



a^oo 48V J a (2vr) 4 



— K y >(p, k){^k^ a - 7 Q /c7 At ) 



d A p ~ 



= _ lim 

a^oo 48V J a (2tt^ 4 



Tr 



p 



lim i6 ° 



a-oo 48V 7a (2tt) 4 



d 4 p ~ 



Tr 



P 4 



K^{p,k)h,k 



(71) 
(72) 
(73) 



with 



10 The authors of [56| do not address the dynamical origin of the singular critical behavior {e.g., 
D^SB), which would be required for the realization of an UV stable fixed point in QED. 



19 



K {2) (p, k) = K {2) (p,p + q,k + q) 

d 



q=0 



K^ a (p,k) EE 
K£ )a (p,k) = 



dq a 
d 2 



R( 2 \p,p + q,k + q) 
K { * 2 \p,p + q,k + q) 



q=0 



dq a dq a q=0 
The trace over spinor and flavor indices is denned as 

tr [L(p)K{p, k)R(k)} = L dc (p)K cd , ab (p, k)R ba (k), 



(74) 
(75) 
(76) 

(77) 



with L and R some projectors and with appropriate summation over double spinor (4) and 
flavor indices (N). 

The BS kernels in Eqs. (]7Tf) — Q73[) contain, in principle, all diagrams except those cor- 
responding to vacuum polarization corrections, since all internal photon propagators are 
canonical or "quenched." As was pointed out in Sec. [TV], the 1/N expansion states that 
the planar diagrams for the a and ir exchanges are dominant. The approximation for the 
BS kernel K^ 2 \ which generates the entire set of planar scalar and pseudoscalar skeleton 
diagrams including ladder photon exchanges for the vacuum polarization is the following: 
the BS kernel is approximated by its "lowest" order skeleton graph, i.e., 



(2) 

K^ d ; ab (p,p + q,k + q) 

kl.ij 



F S cb(k + q,p + q)A s (k - p)^Sad(p, k) 
p)T Pad (p, k) 



(78) 



+ T Pcb (k + q,p + q)A P (k 
+ ^a^kjladlcb D ^{k-p). 

In the symmetric phase, the decomposition of the scalar Yukawa vertex is given by 
Eq. ([Zip. Furthermore, due to chiral symmetry, we have the identities 

(t l5 )T s (k,p) = T P (k,p), A s (q) = Ap(q), (79) 

the a and 7r propagators are degenerate. Second, it was shown in Ref. jlT] that the structure 
function F 2 is rather small compared to the leading structure function Fx (it is assumed that 
Fi describes the leading asymptotic behavior of the Yukawa vertices). Therefore, we neglect 
contributions related to the scalar structure function F 2 . Although it might be possible that 
the contribution coming from gauge interactions is smaller than corrections resulting from 
this structure function F 2 , we keep the gauge interaction in order to compare with results 
mentioned in the literature. Thus, we take for 

5a5k 



kl.ij 



^-^Fx(k + q,p + q)F 1 {p, k)A s (k - p) [l ad l cb + hsadhbcbl 
e o 



+ 8 u S kj D ia/ {k-p)i^rf l 



(80) 

With this truncation for the BS kernel K^ 2 \ we can actually compute the (f)j functions 
(ffl|)- (|7^) , and subsequently analyze the (3 function (|70|) . The truncation fl80D generates an 
infinite series of planar contributions to the vacuum polarization as the leading order in 
1/N, see Fig. |7]. As was discussed in Sec. |IV B| and shown in Appendix [S] up to two loops, 
the scalars and pseudoscalars give the same contribution in the functions <pj. The trace 
over flavor indices yields an overall factor of N in the expressions for cf)j for contributions 
corresponding to a and 7r exchanges. 
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VII. COMPUTATION OF THE (3 FUNCTION 



In computing the functions 0i, 02, an d 03, we initially neglect the ladder photon exchange 
given in Eq. fl80|). Since such contributions were already computed in Ref. |p6| , it will be 
rather easy to include them later in the analysis. 

It is straightforward to show that 02 vanishes, within the proposed approximations. 
Using Eq. (|30D, we obtain from Eqs. (|72|) and (75) that 



Tr 



p^pK {2)a (p, k)(ju,kj a - 7 a fc7p 



oc Tr 



Pl^PilJ^la ~ 7«^7m) 



0. 



(81) 



Thus 2 (ao,0c) = 0. 

With Eq. (0), Eq. (^TJ) for fa reads 



h{oc , 9c) 



— lim 

A^o 

x/tt 



N N 



C 48iV U <W5£ W.CMf.ta *)A»(t " P) 



( 7 ^ 7 Q - 7 a p7 M )(7 M A;7 a - 7aHJ 

+Tr [(j^py" - ^fn^ilsilj^la ~ 7c«^7m)*75 
g? 4 £> (p + k) ■ k 



lim 2Az 

A^oo 



a (2tt) 4 (p + A;) 4 



-[i^+p^A^p), 



(82) 



where we have performed a "harmless" shift of integration, and used the fact that F\ is 
symmetric in the fermion momenta, Fi(p,k) = F±(k,p), because of C-PT invariance. The 
overall factor N results from tracing the flavor Kronecker 8 functions, which is equivalent 
to closing index lines, see Sec. |1V H. After a Wick rotation 



lim — - 



dp 2 



2 f dtt p p 2 (p- k + k 2 



2tt 5 



(p + k) 4 



-[F^k+p^tfAsip). 



(83) 



Since the integrals for the functions <pj are finite, we can write 

oo 

<f>i(a ,g c ) = —N J duG^u), 
o 

where u = p 2 /k 2 , and 



(84) 



A^oo SlT z J 



dQ p k 2 p 2 (k -p + k 2 



2tt 2 



(k + P y 



■[F^k + p, k)] 2 A s (p). 



(85) 



The function G\ is defined with a minus sign to make it a positive function, as will be shown 
to be the case later. The angular integral can be performed if we make use of the following 
Chebyshev expansion: 



The integral is finite, therefore translationally invariant. 
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k 2 (k-p + k 2 ) ^2, 2 2 fc-p 
= c n{k ,p )U n {cosa), cos a - 



(k + P y 



n=0 



kp 



(86) 



where 



7T 



c n (fc ,p ) — — da sia aU n (cosa) 



k 2 (k-p + k 2 ) 
(k + p) 4 



c n (k ,p ) 



(2 + n)0(r-^) f -n6(p 2 -k 2 )[- 



k 



n+2 



(87) 
(88) 



The Chebyshev expansion for the function F\ was already introduced in Eq. 



(and 



Thus, following analogous derivations in the Appendix of Ref. j47|, the function G\ can be 

expressed as 



G x {p 2 lk 2 ^ 



hm^M £ C lmnCl (k 2 ,p 2 )f m (k 2 ,p 2 )Uk 2 ,p 2 ) 

l,rn,n=0 



(89) 



where the fully symmetric index Ci mn — 1, if I +m + n — even and a triangle with sides I, 
m, n exists, i.e., \l — m\ < n < I + m. Otherwise Ci mn = 0. 

We approximate G\ by keeping only the lowest order term in the Chebyshev expansion, 



G^/k 2 " 



pa — lim 



. P 2 A s (p) 



c (k 2 ,p 2 ) f (k 2 ,p 2 ) 



(90) 



where fo is decomposed into the two channel functions Fuv and Fir, see Eq. (j2~4"l). Then, 
the asymptotics, k 2 ^> p 2 , respectively, p 2 ^> k 2 , of G\ are well approximated by the lowest 
order Chebyshev term (|90|). Again, this is the two channel approximation for the Yukawa 
vertices WJ\. However, for momenta k 2 ~ p 2 the channel approximation is not necessarily 



valid. So, how about Gi(l)? Since, from the appendix of Refs. J47|, it follows that the 
Chebyshev coefficients 



/ 2 n(r,/) >o, f 2n+x {k\p l ) <o 



and from Eq. (|88f) that 



c n {k , k ^ 



■l) n 



c 2n (k 2 ,k 2 ) >0, c 2n+1 (k 2 ,k 2 ) <0. 



(91) 



(92) 



ClmnClik ,k )fm(k ,k ) f n (k ,k 

Lm,n=0 



Hence, by taking into account the properties of Ci mn , we conclude that all terms of the series 

(93) 

of Gi are positive, and the lowest order term gives a lower bound on the series, 

2 oo 

c (k 2 ,k 2 )\f (k 2 ,k 2 )} < J2 Q mnCl (k 2 ,k 2 )f m (k 2 ,k 2 )f n (k 2 ,k 2 ). (94) 



l,m,n=0 
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Therefore, the approximation Eq. (|90"D is reliable for the asymptotics k 2 ^> p 2 , and p 2 ^> k 2 . 
Moreover, Eq. (|90|) is a lower bound on Eq. (|89|) at k 2 = p 2 , so that at least we will 
not overestimate the contribution of scalar and pseudoscalar composites to the vacuum 
polarization. 

The function G\ can now be computed, since fo is expressed in terms of the channel 
functions -F\jv and Fm of Eq. ([24]) . Furthermore, from Eq. (|8l) we see that 



c (k 2 ,p 2 ) = 6(k 2 -p 2 ), (95) 



and the only nonzero contribution to G\ of Eq. (|90|) comes from the momenta k > p . 
Thus, using Eqs. (p0|), (|24l) , and (^), we find 

G^A 2 ) « - Jim P ^^[F VY (k 2 ,p 2 )} 2 e(k 2 -p 2 ). (96) 

As was discussed in Sec. |TJ, the ultraviolet channel function F\j\(k 2 ,p 2 ) is proportional to 
(A/p)' ? / 2 and p 2 A$(p) is proportional to (A/p) _r? . Therefore, the cutoff dependence cancels 
in Eq. (|96| ) as was expected and the angular integral Eq. fl85|) can indeed be written in terms 
of a function which depends only on the ratio of p 2 /k 2 . 

The scaling form (p 2 <C A 2 ) for the scalar propagator A s (p) at g = g c (m a = 0) is given 
by Eq. ( p6[) and Eq. (|28|) . The scaling form for the channel function Fjjy(k 2 ,p 2 ) is given by 
Eqs. (|1 and ©, with p 2 < k 2 < A 2 . Inserting Eqs. © (with m a = 0), @, and © 



in Eq. (|9q), we obtain for Gi 

r(2-cu)r(2 + w) 



8u;7(u;)7(— c<j) 

x it 



2 0(l-u), (97) 



where w = p 2 /k 2 . Thus Eq. ( jS4[) is 

M<*o,9c) « -ACiM, (98) 
i 

^(a ) = / rf«Gi(«) > 0. (99) 
o 

The function Gi is positive, hence 0i is negative. The integral over the function Gi can be 
done explicitly by making use of the integral identity 2.15.19.1 in Volume 2 of Prudnikov et 
al. The result is 

n \ = 1 r(i-^) 7 (-^) (x y 

^ l{ao> u 2 \(2+o;) r(l+w) 7 (a;) ^2 J 

x 2 F 3 (2 + u, 1/2 + u; 3 + u, 1 + 2cu, 1 + cu; 2A ) 
- 2 F 3 (2, 1/2; 3, 1 + <j, 1 - u;; 2A ) 

i r(i + ^H (XoY" 

+ (2-u)T{1-u)j{-lo) \ 2j 

X 2 F 3 (2 - w, 1/2 - w; 3 - w, 1 - 2cj, 1 - w; 2A ) 1. (100) 
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The above analysis of the function <pi is repeated for the function 03. The second 
derivative, K^ a (p, k) (see Eq. (|76f)), of the BS kernel given by Eq. (^(]) is 

K^(p, k) oc lim Q^-m + q,p + q). (101) 
The SDE for Fi (in quenched-ladder approximation, see Fig. ||) is given by 

F 1 (k + q,p + q) = 1 - iX Q ( ^ i r2 , + (fe w P) " r) ^ ifr + k-p,r), (102) 

ja 7r r 2 (r + k — p) 2 {r — p — q) 2 

where we recall that we neglect the vertex function F 2 . Thus 

d 2 . f d 4 r (r 2 + (k - p) ■ r) _, . , . 

hm , a aa F 1 (k + q,p + q) = -i\ / -^^7 — —. r^— F x (r + k — p,r) 

q^o oq a oq n Ja it 2 r 2 (r + k — p) 2 



d 2 1 

x lim - - -. (103) 

q^o oq a oq a (r — p — q) 2 



By making use of the identity 

d d 1 



dq a dq a q 2 

we obtain 



4vr 2 ^ 4 (g), (104) 



kl.rj Cq \ p K I 

in Minkowsky formulation. Inserting the above expression in Eq. ([T3|) the equation for 03 
takes the form (in Euclidean formulation) 

*•<<*■*> ^ " fit /*" / $ + p ' fc »l 2A ^ p »- < 106 > 

Then 

00 

<p 3 {a ,g c )=N JduG 3 {u), (107) 


G 3 (P 2 A 2 ) - lim /^ ^•P + f) 3 [ir 1 (fc + p,fc)] 2 A 5 ( P ). (108) 

A^oo 7T 8r J 27T^ (k + p)° 

We use the following Chebyshev expansion: 

(k p + k ) _ y, d n (k 2 ,p 2 )U n (cosa), cos a = ^— ^, (109) 
{k+p) b ^ Q kp 

where 
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d n (k 2 ,p 2 ) 
d n (k 2 ,p 2 ) 



-J da sin aU n (cosa) ^ - fc)6 ; 



(110) 
(111) 



-1) 



n + 1 + 



n + 1 



n-l 



6+E( 4 + 
z=o 

n-l 

E(2-0 



U=0 



3 K -*■>(*; 



n>l. (112) 



The function G3 can be expressed as 
G 3 (P A ) = lim r-^— 



E Q^d^^^Uk^P^Uk 2 ^ 2 ). 



(113) 



i,m,n=0 

We also approximate G3 by keeping only the lowest order term in the Chebyshev expansion, 



G,{p 2 /k 2 ) « lim.-?^#rf (fc 2 , P 2 ) f/ (*V) 



1 2 



(114) 



A^oo 7r 87T" 

Then, again the asymptotics, k 2 3> p 2 , respectively, p 2 ^> k 2 , of G3 are well approximated 
by the lowest order Chebyshev term. Moreover, for momenta k 2 = p 2 the approxima- 
tion Eq. ( |114j ) is exact, since d n (k 2 ,k 2 ) = for all n > 1. Therefore, the approximation 
Eq. ( |114j ) is even better than the analogous approximation, Eq. fl50p, to Gi. Furthermore, 
from Eq. (|111 ) we see that the only nonzero contributions to G3 of Eq. (|114 ) are given by 
momenta k 2 > p 2 . Thus, using Eqs. (|114j) , (24) , and (|111|) , we find 



G 3 (p 2 /k 2 ) « lim (1 - ^^[F vv (k 2 ,p 2 )] 2 9(k 2 -p 2 ). (115) 



Substituting Eqs. (^), (0), and fl3§) in Eq. ( [TToD , we obtain for G 3 



a 



GM = — 1 - — 



7T 



3w\ r(2 - w)r(2 + w) 



4 / 8ury(u;)7(— a;) 



x it 



0(1 -u), 



where u = p 2 /k 2 . Thus Eq. ( [TUT] ) is 



03(ao,fiO ~ NC 3 (a ), 
1 

C 3 (a ) = | duGsH > 0. 




(116) 

(117) 
(118) 



The function 3 is positive, and can be computed in the same way as <f>\. The result is 



(3(^0) = — [Cl(«o) - t(oq)] 
7T 



;ii9) 
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where 



r(a ) 



3A f 1 Tjl-uM-u) (XoV 
4lu 2 \(3 + w) T(1 + lu)^(lu) V 2 ) 
x 2^3 (3 + oj, 1/2 + u; 4 + w, 1 + 2cj, 1 + cj; 2A 

- jj 2 F 3 (3,i/2;4,l + u;,i-cu;2A ) 

i r(i + ^) 7 M /AoV" 
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(120) 



In the computation of the functions 0i, 02, and 03 the ladder (planar) photon exchanges 
have been neglected. After reinstating the ladder photon exchange term of Eq. (|80|), we 
obtain, together with Eqs. (p8|) and (|117|) , that 



0i(ao,^c(ao)) 



2tt 



-iVCi(a ), ^(ao.ffcM) = 0, 03(ao,^c(ao)) = -^Cs(ao)- ( 121 ) 



The ladder photon exchange only contributes to (pi, see again |p6| . After substitution of 
Eq. in Eq. (0), the (3 function reads 



71 



2 

3 + T 



a /2n - Nd(a ) 



N( 3 (a ) 



« /2vr + iVCi(ao) 

where explicit expressions for d and ( 3 are given by Eq. (|100|) and Eqs. (|119|) and (|120|) 



(122) 



VIII. UV STABLE FIXED POINTS 



Let us start analyzing Eq. ( |122j ) by first considering the properties of the functions Ci( a o) 
and C3( a o)- For «o small, the expansion of the functions £i and ( 3 can be computed from 
Eqs. ([T00|) and (|TT9|). The result is 



Ci(«o) ^ + «, CsM « ^ + °( a o), (123) 



showing that £3 vanishes faster that Ci for «o — > 0. The functions Ci( a o) and Cs( a o) have 
been plotted versus a /a c in Fig. [|. First, it is clear that Ci and ( 3 , are positive, and have 
a maximum at some intermediate value of < olq < a c = 71/ 3. For instance, Ci has a 
maximum d « 0.123 at a /a c ~ 0.58 (a; ~ 0.65). Second, the functions Ci and ( 3 vanish at 
the pure NJL point a = in accordance with Eq. ( |123[ ), and at the CPT point a = a c . 
At a — 0, we can consider this is as a reflection of the fact that hyperscaling breaks 
down due to logarithmic corrections; the "effective" Yukawa coupling is trivial, therefore 
vanishes. At «o = oc c , where the critical exponents become singular, the vanishing of £i 
and £3 is related to the dynamics of the conformal phase transition (CPT), which has been 
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thoroughly discussed in Ref. [plj] . There are no light a and ir exchanges in the symmetric 
phase [47 1 which consequently implies the absence of effective Yukawa interactions.^] 

Let us compare the (3 function ( |122| ) with the (3 function (p^) of the gauge-Higgs- 
Yukawa model ( [50D in the 1/N expansion. Then, the entire set of planar a and ir exchanges 
is generated by the kernel 



K { 2 ab {p,p + q,k + q) « (5 ii 5 H ^|-A 5 (A; - p) [l a£ zl c ft + «75 a d*7s 



(124) 



where A^d?) = 1/p 2 . With such a kernel, 2 and 3 are zero, because the right-hand side of 
Eq. ( |124j ) does not depend on the momentum q. The expression for 0i, in this case, can be 
computed straightforwardly (Ay = g y /47r); 



, r OAT 2 • f d P P- k 1 



NX 



Y 



2tt 



(125) 



Again we introduce the ladder photon exchanges by the replacement 

a N\ Y 



0i(Ay) — > 0i(a o ,Ay) 
Hence, in this case, the (3 function is 

f3 a (a , Ay 



2n 



2tt 



Na 2 Q 
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a /2n - NXy/2it 



2 

3 + 1 - a /27r + NX y /2tt 



(126) 



(127) 



Comparing the (3 functions (|122|) and (|127|) leads to the suggestion that Ci( a o) is analogous 
to the Yukawa coupling Ay in a gauge- Higgs- Yukawa model, Ci( a o) ~ Ay /2n. 

This is a crucial point. The general consensus is that for a gauge-Higgs- Yukawa model 
the Yukawa interaction Ay is trivial, thus Ay — > in Eq. (|127|) . However, the situation is 
essentially different for d in the GNJL model. There the "effective" coupling (i is formed by 
the exchange of a and n bosons, with the Yukawa vertices, and (pseudo) scalar propagators 
fully dressed (i.e., the skeleton expansion). The cancellation of the Z factors, see Sec. |TJ, 
which is related to the fact that the hyperscaling equations are satisfied, gives rise to a finite 
nonzero Ci( a o) at the critical curve (go = 9c) for < cto < a c - The other nonzero function 
(3 results from taking into account fully dressed Yukawa vertices. 

Let us now the discuss the possible existence of UV stable fixed points. A necessary but 
not a sufficient condition for the realization of an UV stable fixed point is that N(i has to 
be larger than both N( 3 and ao/27r, and iV£i ~ 0{1). For large N, the contribution of the 
planar photon exchanges (represented by the cko/27t terms) is negligible with respect to iV£i 
and iV£3. Moreover Fig. [| shows, for «o small, that (1 is considerably larger than £3. This 
means that only for flavors iV larger than some critical value N c UV stable fixed points can 
be obtained. 



Moreover at the CPT point four-fermion interaction are marginal instead of relevant, and start 
to mix with the gauge interaction, hence the analysis becomes considerably more complicated. 
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By substituting the expressions ( |100|) and (|119|) for (i and (3 in Eq. (|122|) , we can 
straightforwardly analyze the f3 function graphically. In Fig. |^ the (3 function is plotted for 
various values of N. Figure | shows that for values of iV > N c , with 55 > N c > 54, UV 
stable fixed points exist, the largest being a+ ~ 0.13; 

N = 55 : /3 a (0.13, <? c (0.13)) « 0, v« = -(3' a (0.13, </ c (0.13)) « 0.07, (128) 
N = Q0: j3 o (0.1,j c (0.1))«0, v* = ~P' a (0.1, &(0.1)) « 0.15. (129) 

In accordance with Eqs. fllS]) and (53|), the fixed points are first-order zeros of f3. The general 



pattern is clear; the larger N, with N > N c , the smaller will be the UV stable fixed point, 
but the larger will be the critical exponent rj a . 

The above pattern also suggests that when N — > oo the UV stable fixed point a* — > 
and we would obtain an asymptotically free theory. This is not the case. It was shown in 
Refs. [|45"|,f47| as «o goes to zero that a logarithmic correction appears in the expression for 



the scalar propagator As- In fact, the scaling form for As (with q/A <C 1) is only valid for 
values of u so that (q 2 /A 2 ) w ^> q 2 /A 2 , see Ref. [£7|. The logarithmic correction gives rise to 



the breakdown of hyperscaling relations and is synonymous to triviality of the four-fermion 
interactions (the NJL model). Since our results rely heavily on the existence of scaling forms 
such as Eqs. (|26]) and (0), we can only trust our results for values of a* not too small. 

Since we have made use of results obtained in the quenched approximation, we mention 
that the plots of the (3 function are (at the most) reliable at or in the vicinity of the UV 
stable fixed points at which the quenched approach is self consistent, see Eq. (|6~8"D. 



In Fig. nj the case of iV = 60 fermion flavors is compared with the one-loop (3 function 
of QED. For very small values of «o < 1/100 indeed the one-loop QED result coincides with 
that of the GNJL model, however for larger values of «o the (3 function ( |122| ) deviates from 
the one-loop expression, and eventually an UV stable fixed point is realized at a* ~ 0.1. 

The analysis shows that a rather large number of flavors, 

N > N c « 54, (130) 

is required to obtain UV stable fixed points. From the point of view of the 1/N expansion 
this seems a consistent result, since other than planar contributions are suppressed by at 
least factors of (say) 1/N C . However, from the phenomeno logical point of view, the result is 
unsatisfactory, since it implies that the unquenched Abelian GNJL model (exhibiting UV 
stable fixed points) is only practically applicable for models which have at least N c fermion 
flavors (fractions rounded up). Therefore, it is appropriate to discuss how N c depends on 
the approximation. 

First, we stress that the second term on the right-hand side in Eq. (|122|) containing the 
Ci function causes the suppression of charge screening and is responsible for the possible 
realization of an UV stable fixed point. The denominator in the second term is a direct 
consequence of the resummation of the infinite ladder a and 7r exchanges, and it is mainly 
due to this denominator 1 + N(i that the critical number of fermion flavors N c is large. 

Second, the existence of an UV stable fixed point for a specific number of fermion flavors 
iV depends on the interplay between the functions £i and ^3; which are given in terms of 
integrals of the functions G\ and G3. Let us recall that the lowest order Chebyshev expansion 
for Gi, Eq. (|9~0"D is a lower bound on G\ of Eq. (j§9|), since all terms of the Chebyshev 
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expansion are positive at k 2 = p 2 , the same cannot be said about the approximation (|114| ) 
for G3. Thus keeping more terms in the Chebyshev expansion leads to an increase of (1, 
whereas the effect on ( 3 is less clear, because of the alternating Chebyshev series for £ 3 . 
Therefore, an improvement of the computation of £i will probably lead to a decrease of the 
critical flavor number N c . 

Moreover, in the computation of the functions £i and £3 we have used Yukawa vertices 
(Ts) and a and tx propagators (As) which were obtained in the quenched-ladder approxi- 
mation. An interesting question is whether the improvement of the ladder approximation 
for the gauge interaction (e.g., by including crossed photon exchanges) leads to a increase 
of £1, and thus a decrease of N c . 

Finally, we recall that we have neglected the effect of the Yukawa vertex function F2 
(Eq. pl|)), but clearly the inclusion of F 2 in the analysis could change the results quantita- 
tively. Whether such an improvement will tend to increase or decrease N c remains unclear 
at this stage. 



IX. CONCLUSION 

There are strong indications that four-fermion interactions become relevant near the 
chiral phase transition in GNJL models in four dimensions, due to the appearance of a large 
anomalous dimension. The main objective of this paper was to study the effect of such 
relevant four-fermion interactions on the vacuum polarization of the gauge coupling and to 
reinvestigate the problem of triviality for a particular Abelian GNJL model with N fermion 
flavors. To obtain new results, the four-fermion interactions had to be taken into account 
beyond the commonly used Hartree-Fock or mean-field approach. 

The crucial feature of the GNJL model, within the quenched mean-field approximation, is 
that a nontrivial Yukawa interaction (i.e., an interaction between composite (pseudo)scalars 
and fermions) exists for < «o < ct c . The existence of such a nontrivial Yukawa interaction 
requires the cancellation of renormalization constants of the a and 7r fields in fermion- 
antifermion scattering amplitudes such as the BS kernel K^ 2 \ This is analogous to the 
requirement of hyperscaling (see Sec. [TIT]). If the hyperscaling equations are satisfied, then 
only two of the critical exponents are independent, e.g., r\ and 7, see Eqs. (|^) and (p7|). 
The existence of hyperscaling relations between the critical exponents is intimately connected 
with the existence of Ward-Takahashi identities (and thus the Goldstone mechanism) arising 
from the continuous symmetries of the model. 

The skeleton expansion for the BS kernel provides a natural framework to take into 
account the anomalous dimensions of Yukawa vertices and a and tt propagators. Within 
the skeleton expansion, a and tt exchanges are described in terms of fully dressed Yukawa 
vertices and a and tt propagators. The actual computation of the anomalous dimension, 
and the resolution of the scaling form requires a solution of the SDE's for Yukawa vertices, 
and a and tt bosons. 

In previous work such fully dressed Yukawa vertices and a and tt propagators have 



been analyzed in the quenched-ladder mean-field approximation, see Ref. J47| and references 
therein. To make use of these results consistently, we used the following approximations. 
First, we assumed that the bare coupling parameters are fine-tuned close to the critical 
point, i.e., close to an UV stable fixed point, at which (3 g ~ and f3 a ~ 0. In that case, the 
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quenched or canonical approximation for the photon propagator is self-consistent. Second, 
the gauge-interaction is considered in the ladder form, with bare vertices. Third, we used 
the 1/N expansion (with iV the number of fermion flavors) which states that planar a 
and 7r exchanges describe the leading contribution to Green functions for large N. Then, 
due to the specific form of the chiral symmetry with both scalars and pseudoscalar in the 
adjoint representation of the Ul{N) x Ur(N) symmetry, we argued that in so-called zero-spin 
channels (such as Yukawa vertices and a and 7r propagators) the planar a and 7r exchanges 
cancel each other for momenta larger than the mass of the a boson (in fact in the symmetric 
phase this cancellation is exact). Moreover, an important property of the planar (ladder) 
approximations is that they respect the vector and chiral Ward-Takahashi identities. 

The method of Ref. [fjql , provides a nonperturbative framework independent of the 
fermion wave function Z, and allowed us to compute the contributions of the infinite set of 
planar a and 7r exchanges to the vacuum polarization. The result of the computations is 
that the GNJL model exhibits an UV stable fixed point, (3 a {&*, 9c{ a *)) = 0, for any value 
of N that exceeds some critical value N c (N > N c ). This critical number of flavors turned 
out to be N c pa 54. The larger the number of fermion flavors, the smaller the UV stable 
fixed point a* will be, provided N > N c and a* not too small. Since our results are derived 
on the basis of the existence of hyperscaling laws, we cannot extrapolate our results into 
the region where hyperscaling breaks down due to logarithmic violations, i.e., when a* — > 
(N -> oo). 

From a phenomenological point of view, the large value for N c puts questions to the 
applicability of the GNJL. However, we have given a few arguments in the previous section 
suggesting that N c could be rather sensitive to approximations, and that an improvement 
of the approximations and calculations will probably lead to a smaller value for N c . 

The realization of an UV stable fixed point is motivated by the observation that contri- 
butions of planar a and it exchanges to the vacuum polarization, in an Abelian gauge-Higgs- 
Yukawa model, have identical sign, and tend to reduce screening. In analogy, four-fermion 
interactions describe attractive forces between virtual fermion- ant if ermion pairs in the vac- 
uum polarization. 

The conventional leading term in the vacuum polarization is the one-loop correction 
describing the creation of fermion-antifermion pairs. These virtual pairs can be considered 
as dipoles causing the screening; the vacuum is a medium of the insulator type. Such a 
screening is proportional to the coupling «o and proportional to the number of fermion 
flavors N. However, if a particular fraction of the total amount of fermion-antifermion 
pairs created are correlated by attractive four-fermion interactions, represented by a and 
7r exchanges, then clearly these composite neutral states are not capable of screening. The 
negative term N(i in the f3 function ( 122 ) represents the contributions and the attractive 
nature of four-fermion interactions in the vacuum polarization. 

Within the quenched-ladder mean-field approximation, the critical curve and critical 
exponents are independent of the number of fermion flavors. Within our approximation 
scheme, the mechanism of charge screening clearly is flavor dependent, since the total number 
of virtual fermion-antifermion pairs is proportional to N and the total number of composite 
scalars and pseudoscalars grows as 2N 2 . The larger the number of flavors, the stronger the 
effect of four-fermion interactions. The fixed point appears when the virtual pairs completely 
loose their ability to screen. 
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The existence of an UV stable fixed point implies a nontrivial continuum limit of the 
Abelian GNJL model. The analysis presented here suggest that in the full unquenched 
GNJL model the critical line is replaced by an UV stable fixed point (on the critical line) 
whose exact positions depends on the number of fermion flavors. If the number of fermion 
flavors is below some specific value, the critical four-fermion dynamics are not sufficient to 
yield an UV stable fixed point. In that case the unquenched GNJL model only has a trivial 
(IR) fixed point and the chiral phase transition is of the mean-field type. 
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APPENDIX A: TWO-LOOP VACUUM POLARIZATION 

In this appendix we compute two-loop vacuum polarization corrections including a and 
7r exchanges, see Fig. |5[ We derive the two-loop contribution by making use of the one-loop 
computation of the photon-fermion vertex [|6!^ ,|63[j. 

The SDE for vacuum polarization tensor reads (N = 1) 

IPV) = *eg f Tr [>fS(k + q)r»{k + q, k)S(k)) . (Al) 

Assuming that the WTI's are respected, the vacuum polarization tensor is transverse: 
W v (q) = (-g^q 2 + q^q u ) II(g 2 ), so that 

n(g 2 ) = -^f / A (0 Tr M(* + «)r"(* + ?. k)S(k)\ ■ (A2) 

Let us write and denote the one- loop vertex and self-energy corrections with a subscript (1) 
as follows: 

V(k,p) = 7 M + A{ t 1) (M, S{p) = ^ [l + Z (1) (p 2 )} . (A3) 

Besides a photon exchange, we take into account a scalar and pseudoscalar exchange in the 
one-loop vertex, and self-energy, i.e., 

A^(k,p) = A» lv) (k,p)+A» ls) (k,p)+A» lp) (k,p), (A4) 
pZ (1) (p 2 ) = -E (1V )(p) - E(is)(p) - S(ip)(p). (A5) 

With one-loop vertex corrections 
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—-^(-ie )j x iS(k-w) 



-ie )A^ ls) (k,p) 



A (2vr) 4 

x(-ie )^iS(p - w)(-ie )-f a iDxa(w), 
' —ig Y )liS(k — w) 



-ie )A^ 1P) (k,p) 



a (2tt) 4 

x (-ie )^iS(p - w)(-ig Y )liA s (w) 
igyjilniSlk - w) 



and the self-energies 
«£(iv)(p) 

iE(iP)(p) 



d 4 k 



a (2tt) 4 

x (-ie )jHS (p - w) {-ig Y )ilhi A P (w) 



-seoh^AOHeob^^ " p), 



a (2tt) 4 
d 4 k 

fttxx (-igY)liS(h)(-ig Y )liA s (k - p), 

A 

d A k 



(2tt)< 



-igY)ij5iS(k)(-ig Y )ij5iAp(k - p). 



(A6) 
(A7) 
(A8) 

(A9) 
(AlO) 
(All) 



Taking free massless fermion, scalar, and pseudoscalar propagators, and the photon propa- 
gator in the Feynman gauge (a = 1), 



S(p) = 4, D, v {q) = -ilf, A s (q) = A P (q) 

the one-loop vertices can be expressed as 

Aft v) (M = 2e 2 R»(k,p) + 2elS»(k,p), 
Al ls) (k,p) = -&Rr(k,p)+&S»(k,p), 
^(iP)(k,p) = A^ is) (k,p), 

where the last identity is obtained form Eq. ( |A8| ) by using 757^ = 



1 



(A12) 



(A13) 
(A14) 
(A15) 



-7^75, and where 



R»(k,p) = -i [ 

J A 

S»(k,p) = -i I 

J A 



d 4 w 7 M (p — w)(k — w)/2 — [k — w)(p — w)-y >1 /2 



a (27r) 4 (k — w) 2 (p — w) 2 w 2 

d A w [ (k — w) ■ (p — w)-y^ {k — wY{p — w) 
a (27r) 4 \_{k — w) 2 {p — w) 2 w 2 {k — w) 2 {p — w) 2 w 2 
{p — wY{k — w) 



(A16) 



(k — w) 2 (p — w) 2 w 2 . 
Thus the sum of one-loop vertex corrections, Eq. (A|), can be rewritten as 



Af 1) (k,p) = 2 



2 2 
e o - 9y 



R»(k,p) + 2 e 2 + g 2 Y S"{k,p). 



(A17) 



(A18) 



The sum of self-energy contributions, Eq. 



, can be computed straightforwardly 
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z m (p 2 ) 



16tt 2 



In 



— p 2 



+ 



(A19) 



in Minkowskian formulation. 

The vacuum polarization up to two-loop corrections can be expressed as 



II(g 2 



where 



n ( i)(g 2 ) 



n (2 a)(g 2 ) 



n (2fo) (g 2 



Kir 



n (1) (g 2 ) + n (2a) ( g 2 ) + n (2fe) ( 



#k Tr P^{q)Y{k + q)^k 



3q 2 J A (2tt) 



(k + q) 2 k 2 



K7 



d±k Tr lfi (k + q)^k 



3q 2 J a (2tt) 



(k + q) 2 k 2 



Z (1 Mk + q) 2 ) + Z (1) (k 



d *k Tr^k + ^A^k + q,^ 



3q 2 J A (2tt) 



(k + q) 2 k 2 



(A20) 

(A21) 
(A22) 
(A23) 



The one-loop vacuum polarization Ilm can be computed straightforwardly by making use 
of the projector P^(q) = 9 V — 4%?k I l( t ! > which by contraction with the vacuum polarization 
tensor eliminates the term in proportional to the g^ v tensor. With this projector the 
quadratically divergent contribution to IR^,, which is an artifact of a hard-cutoff regulariza- 
tion, is eliminated explicitly.^ The result is the well-known one-loop vacuum polarization: 



n ( i)(g 2 ) 



«0 

3tt 



In 



A; 

a 2 



+ 0(1) 



(A24) 



with q 2 the Euclidean momentum and chq = ^q/Att. 

The sum of the one-loop vertex functions is given in terms of the functions W 1 and S^. 



The one-loop vertex in the Feynman gauge has been computed in Ref. |62], see also Ref. |63 
for arbitrary gauge. One can show that, with q 2 = (k — p) 2 , 



2e 2 R^k,p) = A» 1R) (k,p), 
2elS^k,p) = Al 1L) (k,p) + Al u) (k,p), 



with 



A?U)(M 



r 1 r 



Z {1) {k 2 



(k + pY(k + p) 
2(k 2 -p 2 ) 



Z {1) (k 2 



z ( Ap 2 ) 



2 

r 2 (A;V,g 2 )TAM + T3^ 
T 8 (k 2 : p 2 ,q 2 )Tg(k,p), 



(A25) 
(A26) 



(A27) 

(A28) 
(A29) 



13 The quadratically divergent contribution A 2 /q 2 is a notorious artifact of computing vacuum 
polarization corrections in the presence of a hard cutoff (i.e., an explicit cutoff in the momentum 
integrations instead of Pauli-Villars regularization, see for a recent discussion Ref. |p4| ). 
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where is the one-loop longitudinal part of the vertex, and where the t{F{ s are the 

one-loop transverse parts as defined and computed in Refs. [|S^ , |63]| . By construction, this 
Ball-Chiu expression for the longitudinal vertex A^^ satisfies the WTI ||62|| : 

= -kZ {1) (k 2 ) + pZ (1) (p 2 ). 
Using Eqs. (|A18|) , ( TA23D , ( |A25| ), and (|A26|) , we write 



(A30) 



n (2 6)(g 2 




where 



n ( 2,)(g 2 ) = - 



n (2L) (g 2 ) + n (27) (g 



d * k Tr^k + fiA^k + q,^ 




{2R) 



(A31) 



3q 2 J A (2tt) 



L, J, 



(A32) 



Since the one-loop transverse vertex functions themselves are finite, i.e., these function 
are independent of the cutoff A, the leading logarithmic contributions to the vacuum polar- 
ization result from integrations over momenta k 2 3> q 2 in H(2R) and 11(2/). These leading 
logarithmic contribution can be found by first deriving the k 2 3> q 2 asymptotic behavior of 
the transverse structure functions given in |k|, after which the integration over angles can 
be performed. In the Feynman gauge a — 1, the asymptotic behavior k 2 ^> q 2 of the r's is 



T-2 



«o 1 
24vr k v 



T 3 



6tt k 2 



(2k ■ q + q 2 ) a 1 
2 24vr£; 4 ' 




(A33) 



where k 2 and g 2 are Minskowskian momenta. By making use of these asymptotic expressions 
the integration over angles in Tl(2R) and R(2i) can be performed straightforwardly, after 
performing a Wick rotation. The integrations over momenta k 2 > q 2 leads to logarithmic 
corrections. The result reads (in Euclidean formulation): 



n 



(2/) 



(Q 2 ) = ^ 



rWg 2 ) 



7T Z 



On 



7T- 1 



— In + 0(1) 




(A34) 
(A35) 



The logarithmic corrections of r 2 and tq cancel each other, and the contributions of T3 give 
rise to a In 2 term. 

An analogous computation can be performed for the self-energy and longitudinal vertex 
corrections. Due to the Ball-Chiu expression ( |A27[ ) for Aj^, the contributions n( 2o ) and 
ri( 2 L) depend on the one-loop computation of the self-energy Z^ given in Eq. (|A19|) . After 
expanding Zm((k + q) 2 ) for k 2 3> q 2 , 



Z (1) ((k + q) 2 ) « Z (1) (k 2 ) + {2k-q + q 2 )Z' {1) (k 2 ) 

+ ~(2A: • q + q 2 fZ'L{k 2 ) + \{2k ■ q + q 2 fZ'!'(k 



6 



(A36) 
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and using that Z^((k + q) 2 ) ~ Z^(q 2 ) for q 2 ^> k 2 , the angular integration can be per- 
formed, and the logarithmic corrections can be computed. The result is 



n (2a) (g 2 ) + (i + ^)n (2L) ( g 

a 



a (a + A 5 



7T 2 








--In 2 


w) 


24 





144 \ q 2 



0(1) 



(A37) 



with ao = Cq/Att and Ay = g Y /^- Thus, comparing this expression with the 11(27) term in 
Eq. ( |A31| ), we see that the "overlapping divergencies" (i.e., the In 2 ) cancel 



IT (2a) (g 2 )+ i + — 



Ay 



n (2L) (g 2 ) + n (27) (g 



a (a + Ay) 



7T^ 



0(1). 



(A38) 



Such a cancellation occurs in a similar manner in any covariant gauge a. Thus, the two- 
loop contribution contribution to II is described solely by the part of the transverse vertex 
containing the Tg 1 tensor]^] i.e., H(2R), and, after adding all the pieces, we deduce that 



vy ' 2vr V 3 2vr 2tt 



(1 



(a /7r)O(l). 



(A39) 



APPENDIX B: THE JOHNSON WILLEY BAKER EQUATION 

In this appendix we derive the equation 



f(i) 



Na 
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0! + 2 (2 + 2 ) 
1-c/h 



+ 



(Bl) 



for the /(i) function given by Eq. (|64"D, with the functions 4>j given by Eqs. (|7T]) ([73D . The 
derivation of Eq. ( Bl|) was given (for pure QED) by Johnson, Willey, and Baker (JWB) 
in Ref. [59. Since their result is formulated in terms of the BS fermion-fermion scattering 
kernel their method is also applicable to the GNJL model. 
In order to derive the result of pEL the following is assumed. 



The fermion wave function equals one, Z = l/A — 1, in the Landau gauge. In 
principle, this assumption is redundant since the JWB result is valid in any gauge. 

Internal photon propagators are replaced by their canonical form A(g) = 1/q 2 which is 
self-consistent in the neighborhood of an UV stable fixed point. Only a single fermion 
loop, thus a single power of In A contributes to the vacuum polarization. 



14 As was shown in [p3|, this particular transverse structure function rs does not depend on the 
gauge parameter a. 
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Translational invariance of naively logarithmically divergent and finite momentum 
space integrals is assumed. Implicitly, use is made of invariance under charge- 
conjugation (C) and parity-time (PT) transformations. 

Also it is assumed that we are in the scaling region of the theory, where the only 
relevant dimensionless variable is q 2 j A 2 . We consider short distances with respect to 
the IR length scale £ ~ l/|m CT |, thus \m a \ 2 < q 2 <C A 2 . 



Hence, 



Z{k 2 ) = l/A{k 2 ) = 1, S(k) = A r»(k, k) = -f. (B2) 



The vacuum polarization tensor is 

W(q) = jf d'k Tr [S(k + q)T^k + q, k)S(k)Y] , (B3) 

where Tr denotes the sum over both spinor and flavor indices. Since, in the chiral symmetric 
phase, the N fermions are degenerate, all N fermion propagators and N photon-fermion 
vertices are degenerate. Hence, the sum over flavor indices gives rise to a factor N, i.e., 
Tr — > iVTr, where Tr is the sum over spinor indices. 

Since vacuum polarization tensor is transverse and the only relevant momentum variable 
is q 2 /A 2 <C 1, the equation for the vacuum polarization can be written as 

n(g 2 ) = - l /f^^) + °M + ° ((^/ A ) ff ) > ( B4 ) 

where a is some positive power. After inserting Eq. ( p3|) , and setting q 2 = in the integrand, 
and using q as the infrared cutoff in the momentum integral, we obtain 

n(g 2 ) - -^^/^^TrK^r-W^T^ 

+2S a {k)T^' a (k)S(k)Y + S{k)r% a {k)S{k)Y], (B5) 
where the derivatives are defined as follows (with T^^k) = T At (/c, k)): 

S a (k) = ^S(k) = SZ{k) = g^S(k) = -g, (B6) 



d 

T nA k ) = -Qq^ T ^ k + ^ k ^> 



1 'a 1 



q=0 



(B7) 



Since the integral Eq. (|BB|) can only be proportional to g^ u , it reduces to 



w -TZTT I d*kTr[S:(k)T»(k)S(kh, 

VOTY 6 Jq,A 



+2S a (k)T^ a (k)S(k) 1 , + S(k)T% a (k)S(k) lfl ]. (B8) 
The SDE for the vertex T M reads, in terms of (see Fig. 
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5 tJ TUk + q,k) = dy7& + ie 2 J [S(p + q)5 nm T^p + q,p)S(p)] dc 



Ia (2tt) 4 
x K% ab (p,p + q,k + q). 

mn,ij 



(B9) 
with 



From now on we omit spinor and flavor indices. Differentiating now the SDE Eq. 
respect to q, and setting q = for the integrand and q as IR cutoff, we obtain 

iW*) = fcS / 7^4 ^(p)rM,a(p)^(p)^ (2) (p, fc) 

+ s«(p)r M (p)£(p)#W(p, fc) + s(p)r^p)S(p)K^( P , *o] , (bio) 

and for the second derivative of the vertex, we find 

rr(*) = fcg / t^ti kp)rnp)S(p)K {2) b, fc) + 2S Q (p)r^(p)£(p)i^ 2 )(p, fc) 
+^(p)r"(p)s r (p) J K'( 2 )(p > fc) + 2S(p)r»> a (p)s(p)K®(p, k) 

+2S a (p)T^p)S(p)K^(p, k) + S(p)T^p)S(p)K^ a (p, fc)J , (Bll) 

where K^ 2 '(j>, k) and the derivatives Kf®(p, k) and K^ a (p, k) are defined in Eqs. Q74])-(f76D . 

The first derivative of the vertex is antisymmetric in a and /x, because of the assumption 
Eq. (|B2|) . Furthermore, C and PT invariance imply that the only nonzero contribution to 
the first derivative of T 1 ' (with q = and A — ► oo) must be proportional to the tensor 
{ifjkla ~ lakjfj). Thus we write 



A: 2 



r( M , Q )(A;) = IV ia (fc) + V a ^(k) = S^ a (k) 



(B12) 
(B13) 



where T' is a dimensionless scalar function.^ Since S u ^(k) = 0, due to the WTI for the 
vertex and Eq. ([B2|), we find that 



(B14) 



After some algebra (taking traces over spinor and flavor indices), we can derive from 
Eqs. (IBTOD and (gig ) that 



lim r MjQ (/c) 
lim r' 

5^0,A-^oo 



>i + 02 + 0i Am T' 

q— >0,A^oo 



2k 2 



1 - 



(B15) 



15 The function R is related to the transverse structure function Ts(k 2 ,k 2 , 0) of Ref. [63]. At the 
one-loop level R^ = -k 2 T 8 (k 2 ,k 2 ,0). 
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where the 4>j functions are defined in Eqs. ([Tl|)-([73|). 

Since r M (fc) = 7^, and using Eq. ( Bll|) for T^ a (k), the second derivative of T^ a in 
Eq. ( [B8| ) can be eliminated. The result is 



n(g 2 



n=l 



where 



h 
h 
h 
h 
h 



d 4 k~ 



48 Jq,A 7T 2 



Tr [S:(k)YS(k) % 



d 4 k~ 



24 Jg,A 7T 

d^k 



2 TT[S a (k)T^(k)S(k)^}, 



24 y g ,A vr 2 J q ,A (2tt) 



<i 4 p ~ 



Tr 



d^k 



24 y g ,A vr 2 7 g ,A (2tt) 



<i 4 p ~ 



Tr 



d 4 £; 



48 J q ,A vr 2 y g , A (2tt) 4 



d A P ^ 



Tr 



5(p)r M , a (p)5(p)^ 2 ) Q (p,A;)5(fc)7^(fc) 



5 Q (p)7^b)^ (2)a (P^)^(^)7 At 5(A ; ) 
Siph.S^K^^^Sik^Sik) 



(B16) 



(B17) 
(B18) 
(B19) 
(B20) 
(B21) 



Using translational invariance, C-PT invariance, Eqs. (|7T1)-(73) and (|B15|) , we can derive 
that 



h 



A 2 

2 r dk 2 



+ 



1 - 



f dk 2 /n/ . 



+ 



/|r + h = <hf^ + 0(l), J 6 = fc/^ + 0(l), (B22) 



where g 2 is the Euclidean momentum. Substituting Eqs. ( |B22|) in Eq. ( |B16|) , we get 



i% 2 



iVa 
~2tT 



2 1 + 2 ( 2 + 2 ) ' 

q + 1 ^ + ^ 3 

3 l-0i 



\n^ + (Na /n)O(l). 



(B23) 



With Eq. ([64]), we obtain Eq. (|B1|). This is the main result of Ref. [[56]. The entire derivation 
did not yet specify the BS kernel K^ 2 \ Therefore Eq. 
as well. 



is applicable to the GNJL model 
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FIGURES 




FIG. 1. The vacuum polarization tensor H /ll/ (q) = (—q 2 g^ u + <2yt<2v)II((/ 2 ), with the blobs repre- 
senting full fermion propagators and a full photon-fermion vertex. 




FIG. 2. SDE for the scalar propagator or a boson As(q). 




FIG. 3. SDE for the scalar vertex or Yukawa vertex Ts(p + q,p) in the quenched-ladder ap- 
proximation. 
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FIG. 4. Skeleton expansion for the BS kernel K^ 2 \ 




FIG. 5. One and two-loop vacuum polarization corrections containing both an internal photon 
(wavy) and a a and tt exchange (dashed). The contribution of it equals the contribution of a. 




FIG. 6. SDE for the full photon-fermion vertex r^(k,p). 
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FIG. 8. The functions £1 and (3 plotted versus a>o/a c . 
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FIG. 9. Plot of the function (3 a versus the gauge coupling ao for various values of the fermion 
flavor number N. 
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FIG. 10. The one-loop f3 function for N = 60 compared with the (3 function including 
four- fermion interactions. 
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